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MINIFICATION  PROCESSES 

Peter  A.  W.  Lewis 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93940 

and 

Ed  McKenzie 

Department  of  Mathematics 
University  of  Strathclyde 
Glasgow,  Scotlamd 

Abstract 

It  is  shown  that  the  autoregrcMive,  Markovian  minification  processes 
introduced  by  Tavares  and  Sim  can  be  extended  to  marginal  distributions 
other  than  the  exponential  and  WeibuU  distributions.  Necessary  and  suffi¬ 
cient  conditions  on  the  hasard  rate  of  the  marginal  distributions  are  given 
for  a  minification  process  to  exist.  Results  are  given  for  the  derivation  of 
the  autocorrelation  function;  these  correct  the  expression  for  the  WeibuU 
given  by  Sim.  Monotonic  transformations  of  the  minification  processes  are 
also  discussed  and  generate  a  whole  new  class  of  autoregressive  processes 
with  fixed  marginal  distributions.  Processes  generated  by  a  maximum  op¬ 
eration  are  also  introduced  and  a  comparison  of  three  different  Markovian 
processes  with  uniform  marginal  distributions  are  given. 

Keywords:  Minification;  Time  series 

1  Introduction 

In  a  series  of  papers  Tavare8(1977, 1980a, b)  introduced  two  stationary  Markov 
processes  with  similar  structural  form  which  he  had  found  useful  in  hydrological 
applications.  The  first  process,  a  maximum  process,  has  an  extreme  value  or 
Gumbel’s  marginal  distribution,  for  which  the  density  function  is 

fx{x)  =  ^  >  0;  -oo  <  I  <  oo  (1) 
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and  the  hazard  rate  is 


Ax(*)  =  fx{x)/Sx{x)  =  (l  -  e-*-*")  (2) 

where  Sx{x)  =  1  -  fx(®)  =  >  x)  is  the  survivcsr  function.  The  second 

process,  a  minimum  or  minification  process,  has  negative  exponential  marginal 
distribution  for  which 

fx{x)  =  pe-f>‘  x>0J>0  (3) 

and 

Ax(a:)  =  0-  (4) 

This  second  process  is  also  investigated  by  Daley,  Chernick  and  Littlejohn(1988) 
who  demonstrate  some  very  interesting  theoretical  properties,  in  particular  that 
the  process  is- a  time  reversed  version  of  the  linear,  additive  EAR(l)  model  of 
Gaver  and  Lewis(l980).  Further,  Sim(198S}  has  shown  that  the  structural  form 
of  this  minification  process  of  Tavares  will  also  accommodate  a  Weibull(K,^) 
marginal  distribution.  The  Weibull  {k,0)  density  function  is 

fxix)  =  X  >  0;k  >  0;^  >  0  (5) 

and 

Xx{x)  =  kI3'^x'^~^.  (6) 

Thus  the  hazard  rate  is  a  power  law  function,  decreasing  from  infinity  to  zero 
&a  X  increases  when  k  <  1,  and  increasing  monotonically  from  zero  to  infinity  as 
X  increases  when  k  >  1.  For  k  =  1  this  is  an  exponential  distribution. 

This  second  process  is  called  a  minification  process  because  the  observations 
{Xn},  where  the  X„’s  are  positive  valued,  are  generated  by  the  equation 

Xn  =  Kxmn(X„,Zn-i).  (7) 

Here  A"  >  1  is  a  constant,  and  {Z„}  is  an  innovation  process  of  independent 
and  identically  distributed  random  variables  chosen  to  ensure  that  {A„}  is  a 
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stationary  Markov  process  with  marginal  distribution  function  Fx{x).  In  the 
negative  exponential  case  (  3)  it  is  found  that  the  autocorrelation  function  of 
the  process  has  the  familiar  geometrically  decreasing  form  of  an  autoregressive 
process  of  order  one,  AR(1),  i.e.,  pxO)  =  J  =  I»2, _ 

The  purpose  of  this  work  is  to  explore  the  generality  of  {  7)  in  two  distinct 
ways.  First  we  attempt  to  determine  the  range  of  possible  marginal  distribution 
functions  Fx{x)  for  non-negative  Xn-  We  shall  see  that  it  is  possible  to  specify 
fairly  simple  conditions  in  terms  of  the  hazard  function  of  the  distribution  which 
determine  whether  that  distribution  may  be  used  as  a  marginal  distribution  in 
(7).  The  second  approach  to  this  investigation  arises  from  the  fact  that  the 
structure  of  (  7)  is  so  simple  that  many  of  the  important  features  of  the  process 
may  be  invariant  under  instantaneous  monotone  transformation.  In  this  way, 
we  may  derive  simple  Markov  processes  based  on  (  7)  but  with  any  marginal 
distribution  we  wish,  even  if  it  is  one  for  which  a  minimization  process  does  not 
exist.  Moreover  Me  may  also  immediately  deduce  many  of  the  beuic  properties 
of  the  resulting  process  directly  from  our  results. 

Sequences  of  non-negative  random  variables  find  applications  in  many  fields. 
The  work  of  Tavares  was  motivated  by  hydrological  considerations,  for  example 
modelling  of  run-off  data.  This  data  tends  to  have  long  tails  and  thus  cannot  be 
modelled  by  exponential  processes  such  as  the  random  linear  coefficient  EAR(l) 
processes  of  Gaver  and  Lewis  (1980).  (See  Lewis,  1985,  for  a  summary  of  these 
models).  Weibull  or  extreme- value  random  variables  are  commonly  used  for 
modelling  the  marginal  distribution  functions  of  run-off  series,  but  processes  with 
these  marginal  distributions  cannot  be  generated  with  linear  random  coefficient 
models.  Thus  the  minimization  processs  are  important  as  a  source  of  time  series 
for  such  processes. 

Another  case  of  interest  is  time  series  of  wind  velocity  magnitudes.  These 
again  are  positive  valued  random  variables  and  their  simulation  is  important 
for  driving,  for  example,  models  of  temperature  mixing  in  the  ocean.  Brown, 
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Katz,  and  Murphy  (1984)  note  that  although  studies  have  shown  that  Weibull 
marginal  distribution  have  been  found  adequate  for  wind  velocity  magnitudes, 
unfortunately,  “no  time  series  models  have  been  rigorously  developed  for  random 
variables  possessing  a  Weibull  distribution  They  therefore  resort  to  transfor¬ 
mations  of  the  data.  Wind  power  data,  being  the  square  of  wind  velocity  data, 
is  even  more  likely  to  need  very  long  tailed  marginal  distributions. 

Again,  in  reliability  studies,  sequences  of  times-between-failures  are  corre¬ 
lated  and  models  are  chosen  on  the  basis  of  a  generally  non-constant  marginal 
hazard  rate.  We  show  how  it  is  possible  to  generate  such  sequences  with  mini- 
fication  processes,  for  example  with  the  familiar  ‘bath-tub’  hazard  rate.  An 
analytical  representation  for  this  type  of  hazard  rate  is  given  in  Gaver  and 
Acar(1979). 

2  The  General  Minification  Process. 


For  the  moment  we  make  no  assumptions  about  the  marginal  distribution  of 
{X„}  or  of  {Zn}  save  only  to  assume  that  a  distribution  can  be  found  for  Z„ 
so  that  is  a  stationary  Markov  process  given  by  (  7)  .  Suppose  now  that 
the  survivor  function  of  the  non-negative  valued  random  variable  X„  is  given  by 
Sx(i)  =  P{X  >  i).  It  is  easily  verified  from  (  7)  that  the  survivor  function  of 
Zn  must  satisfy 


5z(x)  =  Sx{Kx)/Sx{x)  x>0;  K>1  (8) 


This  shows  that  K  must  be  greater  than  one;  otherwise,  since  Sx{x)  is  generally 
decreasing,  the  function  ^^(x)  would  be  greater  than  one  in  value  for  some  i. 

Note  that  for  the  general  process  we  may  write  the  survivor  function  Sx(®)  in 
terms  of  the  cumulative  hazard  Ax(x),  or  the  hazard  rate,  Ax(x),  thus:  Sx(x)  = 
expf-Ax(x)]  =  exp[-  /g*  AxCO*^*!-  Thus,  equation  (  8)  may  be  recast  in  terms 
of  hazard  functions  as  follows: 


fx  fKx 

Az(x)  =  ~  J  Ax(0<^* 


X  >  0 


(9) 
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We  now  consider  what  is  the  set  of  possible  marginal  distributions  for  Xn  in 
the  general  miniiication  process  (  7).  Clearly,  a  necessary  and  sufficient  condition 
for  a  distribution  to  be  suitable  for  this  purpose  is  that  the  right  hand  side  of 
(8)  is  a  survivor  function,  or  equivalently  that  the  right  hand  side  of  (9)  is  a 
cumulative  hazard  function.  In  the  latter  case,  we  require  that  -ln5z(i)  = 
S.z[^x)  —  Xx{t)dt  be  a  non-decreasing  function  of  i,  for  all  x  >  0,  and  be 
increasing  for  some  x  >  0.  This  is  equivalent  to 

r  {K\x{K8)-Xx{8)}  da  >0 
Jy 

for  all  z  >  y  >  0,  and  positive  for  some  z  >  y  >  0,  which  reduces  to 
KXxiKa)  >  Ax(«)  for  all  8  >  0,  with  inequality  for  some  s  >  0.  Multiply¬ 
ing  by  8  yields  the  alternative  necessary  and  sufficient  condition  that 

xXx(x)  (10) 


is  a  non-decreasing  function  of  z,  for  all  z  >  0,  and  an  increasing  one  for  some 
z  >  0. 

Now  assume  that  (  10)  is  true.  Then  as  z  — ►  oo  either  (i)  Sz(x)  — »  0  or  (ii) 
Sz(x)  -*  p,  where  p€(0, 1). 

In  case  (i)  Sz(x)  is  a  proper  survivor  function  and  so  the  minification  process 
(  7)  can  be  constructed  with  the  required  marginal  distribution  for  {A^n}  by  using 
an  innovation  process  {Zn}  whose  survivor  function  is  defined  by  (  8). 

In  case  (ii),  Sz(x)  is  not  strictly  a  proper  survivor  function  having,  in  effect, 
an  atom  of  probability  p  located  at  infinity.  Such  a  property  would  seem  to 
rule  out  this  case  from  a  practical  point  view.  However,  it  is  a  property  which 
is  readily  interpreted  in  practice  because  of  the  form  of  (  7).  Simply,  we  may 
rewrite  (  7)  in  the  form 


(11) 


_  J  KXn-i  with  probability  p 

(  K  wm(Xn-i,Z*)  with  probability  (1  -  p) 

and  the  “new”  innovation  r.v.  Z*  is  simply  Zn  conditional  on  Zn  <  oo.  Thus, 

Sz(z)  -p 


Sz’ix)  = 


1  -p 


(12) 
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In  addition,  although  the  form  of  (11)  is  different,  of  the  process  to  be  discussed 
below  can  be  derived  similarly  provided  we  work  in  terms  of  Sz{x)  =  p  +  {l-p) 
Sz'(x).  Note  that  from  (11)  can  see  that  sample  paths  of  the  process  {X„} 
will  tend  to  exhibit  a  “runs  up”  type  of  behavior.  This  type  of  behavior  is 
characteristic  of,  for  example,  river  flow  data,  but  the  geometric  increase  implied 
by  (11)  may  be  too  severe  for  general  use. 

3  Bivariate  Distributions  and  Autocorrelations 

Using  (  7)  and  (  8)  it  is  straightforward  to  show  that  the  bivariate  distribution 
of  any  two  values  in  the  process  has  survivor  function 

Sx„.X„-j{y,  x)  ^  P  (Xn  >  y,Xn-j  >  X) 

^P{Xn-l  >  y/K,Zr,  >  y/K,Xn-i  >  x} 
from  (7).  By  repeated  use  of  (7)  we  get 

Sxn,Xn-iiyy  (max  {x,y/Ki))  Sx{y)/Sx  {y/K>) 

_  (  Sx{x)Sx{y)/Sx  y  <  K’x 

\  Sx(y)  y>  K>x 

Note  also  that  this  is  a  not  an  absolutely  continuous  distribution  because, 
from  (  3) ,  there  will  be  a  non-zero  value  for 

Pj  —  P  (Xn  —  X‘Xn-j) 

(14) 

=  P(Zi  >  Xo,  Z3  >  XXo, K^-^Xo) 
which  may  be  evaluated  as 

Hence,  the  bivariate  distribution  has  probability  p,  on  the  line  Xn  =  K’ Xn-j 
and  probability  (1  -  py)  spread  over  the  region  defined  by  Xn  <  K^Xn-j  with 
survivor  function  given  by  (  13).  It  is  important  to  note  that  since  the  process 
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is  Markov  all  distributional  behavior  is  characterized  by  such  bivariate  distribu¬ 
tions,  especially  the  form  for  contiguous  observations,  i.e.  fc  =  1. 

In  addition,  note  that  the  bivariate  distributions  of  (X„,Xn-i)  and  (X„,X„-^) 
differ  only  in  that  K  in  the  former  becomes  in  the  latter.  Thus,  the  bivari¬ 
ate  distribution  of  {Xn,Xn-j)  and  its  properties  are  easily  derived  from  those 
of  (X„,X„_i)  by  replacing  K  by  .  In  particular,  we  are  interested  in  the 
autocorrelation  function  of  the  process  {Jfn}.  Px(j)  =  corr(X„,X„_y),  j  = 
0,1, —  Thus,  if  px(f)  =  then  pxU)  ~  J  =  —  This  is  a 

useful  property  since  we  can  now  derive  the  autocorrelation  function  for  any  lag 
j  from  alone. 

Tavares  (1980)  claims  to  show  that  in  the  case  where  {X„}  defined  in  (7) 
heis  a  negative  exponential  distribution  the  autocorrelation  function  of  lag  j  is 
given  by  px{j)  —  Pj,  where  py  is  defined  by  (14).  In  this  negative  exponential 
case,  pj  =  [1/ Ky .  Sim  (1986)  uses  exactly  the  same  argument  in  the  case  when 
{Xn}  is  marginally  Weibull  and  again  derives  px(j')  =  ■  These  results  are 

of  particular  interest  since  this  geometric  autocorrelation  function  is  associated 
with  the  well  known  autogressive  process  of  order  one,  the  AR(1). 

Unfortunately,  although  the  autocorrelation  result  is  true  for  the  negative 
exponential  case  (as  may  be  seen  from  Chernick,  Daley,  Littlejohn,  1988),  the 
proof  indicated  by  Tavares  does  not  appear  to  be  valid  and  certainly  does  not 
extend  to  the  Weibull  case,  as  Sim  states.  In  general,  px{j)  /  Pj  although 
equality  holds  in  the  negative  exponential  case.  By  considering  E{Xn  | 
we  may  show  that  for  the  general  stationary  process  defined  by  (7) 

£;(A'„A„-i)-X£:|x 

Using  E  I^K  Sz{z)dz\  =  E{X)  —  mx,  we  may  extend  (  16)  to  obtain  an 
expression  for  the  covariance  thus; 

=  KE 


C{X„,Xn-l) 


{X-mx)j^  5z(z)dr|. 
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Hence,  px(l)  may  be  obtained  from  (  17)  by  dividing  by  Var(Jfn),  and  px{j)  may 
•'  obtained  by  replacing  K  in  px(l)  by  ■  (In  this  context,  recall  that  Sz{z) 
is  a  function  of  K).  The  geometric  autocorrelation  for  the  negative  exponential 
case  claimed  by  Tavares  may  be  readily  shown  from  (  17).  In  passing,  we  note 
that  the  autocorrelation  function  has  this  geometric  form  in  generad  if  and  only 
if  px(l)  =  some  a  >  0.  The  quanitity  p,  defined  by  (14)  is  also  a 

useful  measure  of  dependence,  although  we  do  not  consider  it  in  any  detail  here. 
The  sequence  (py  ;  j  =  0, 1,2, . . .}  is  the  survivor  function  of  the  length,  T  say, 
of  runs  of  the  form  (X,  KX,  K^X, . . . ,  K'^X),  and  so  P(T  >  j)  =  pj. 


4  Examples  of  Minificaton  Processes 


(i)  A  simple  example  of  a  distribution  which  cannot  be  the  marginal  distri¬ 
bution  for  a  minification  process  as  defined  by  (7)  is  given  by  taking  X 
to  be  uniformly  distributed  over  (0,1)  with  probability  0.5  and  uniformly 
distributed  over  (1,6)  with  probability  0.5.  In  this  case 


_  f  1  —  0.5i 
“  \  0.6-  0.1* 


0  <  *  <  1 
1  <  *  <  6 


and  it  is  a  matter  of  straightforward  calculation  to  show  that  if  1  <  if  <  3 
then  Sz  as  defined  by  (8)  is  increasing  for  all  z  in  the  interval  (l/if,  l). 
We  may  also  verify  that  the  condition  given  by  (10)  is  violated. 


(ii)  .A.n  important  but  simple  example  of  a  minification  process  is  provided  by 
the  uniform  distribution  on  (0,1).  In  this  case,  Sx(x)  =  1  -  z,0  <  z  <  1, 
and  zAx(3;)  =  z/(l  -  z)  is  clearly  increasing  so  that  (10)  is  satisfied. 
Further,  from  (8) 


1  2^  Xa 

which  is  a  proper  survivor  function  and  Z  is  given  by  Z  -=1] /(K  ~  l  +  U), 
where  U  is  uniform  on  (0,1).  Ip  addition,  the  autocorrelation  function 
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of  {X„}  can  be  derived  from  (ll)  and  is  found  to  be  Pjf(l)  =  1/if,  so 
that  pxU)  =  (l/if)^,  J  ==  0, 1, . . ..  Thus,  this  uniform  process  enjoys  the 
geometrically  decaying  autocorrelation  of  the  AR(l)  process.  Processes 
which  are  marginally  uniform  are  important  since  the  uniform  random 
variable  X„  can  be  given  any  other  distribution  by  means  of  the  inverse 
distribution  function  transformation.  Thus,  if  we  wish  a  random  variable 
Y  with  distribution  functions  ^^(y),  we  use  Y  =  Fy  ^(X).  This  idea  will 
be  discussed  in  more  detail  later. 

(iii)  It  is  clear  from  (10)  that  Xn  defined  by  (7)  may  have  any  marginal  distri¬ 
bution  whose  hazard  rate  is  itself  non-decreasing,  e.g.  the  uniform  above,  a 
Gamma  (k,0)  with  /c  >  1  or  a  Weibull  {k,0)  with  k  >  1.  We  consider  now 
a  distribution  whose  hazard  decreases  over  the  sample  space,  the  Weibull 
(k,/9)  with  /c  <  1.  This  is  detailed  in  Sim  (1986),  but  we  note  here  that 
x\{x)  =  K/J*!*  is  increasing  for  all  <c,  and  so  condition  (10)  is  satisfied. 
Further,  5^(1)  =  exp{-(if*  -  l)(^i)*}  -+  0  as  z  -♦  oo.  Thus,  the  gen¬ 
eral  minification  process  (7)  accommodates  the  Weibull  distribution  for  all 
K  >  0.  However,  as  noted  above,  Sim’s  derivation  of  the  autocorrelation  is 
wrong.  He  shows,  correctly  in  this  case,  that  pj  =  (l/X)^,  but  pxU)  ^  Pj 
in  general.  Using  (17),  it  is  possible  to  show  that  the  form  of  px(l)  for 
genera!  K  cannot  easily  be  derived.  The  case  when  k  =  2  is  tractable, 
however,  and  in  this  case 

'  (i)  +  - 1}/  ('  “  i) 

(iv)  The  Pareto  distribution  provides  an  example  in  which  the  form  (11)  is 
required  rather  than  (7).  Here  Ax(x)  =  “/(I  +  x)  with  a  >  0  a  shape  pa¬ 
rameter,  and  condition  (10)  is  again  satisfied.  However,  Sz{x)  =  {(1  +  x) 
/  (1  +  ifi)}“  —»  K~°‘  as  I  — »  oo.  Hence,  the  process  {X„}  with  this 
Pareto  marginal  distribution  may  be  generated  using  (11)  and  (12)  with 
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p  =  K  ^  and  an  innovation  process  [Zn}  whose  survivor  function  is 

{[Kil  +  x)/il  +  Kx)r-l}/{K‘‘-l). 

Again  the  autocorrelation  function  can  be  derived  from  (17)  and  we  find 
that,  for  a  >  2,px{j)  =  {l/Ky,j  =  0,1, —  Thus,  the  Pareto,  like  the 
uniform  minification  process,  has  the  familiar  geometric  autocorrelations 
of  the  AR(1)  process. 

(v)  An  interesting  case  is  that  of  the  so-called  bathtub  hazard  rate  (Gaver 
and  Acar,  1979),  which  we  could  model  as  A(*)  =  paF x’'’~^  +  ^  +  P'^’ x^~^ 
with  0  <  p  <  1  <  P.  This  models  a  situation  where  components  have  high 
likelihood  of  early,  infant  failure,  otherwise  have  a  constant  hazard  rate 
and  then  finally  reach  a  “wear-out”  state  corresponding  to  the  Weibull 
distribution  with  P  >  1.  This  is  actually  the  hazard  rate  of  a  random 
variable  which  is  generated  as  the  minimum  of  three  independent  ran¬ 
dom  variables,  one  being  exponential(/9),  the  others  being  Weibull(p,a) 
and  Weibul^P,'/).  Then  xA(i)  =  poPx^  +  +  P'/’r^,  which  is  clearly 

increasing  in  x,  so  that  a  minification  process  exists.  In  fact 

Xz{x)  =  KXxiKx)  -  Ax(x)  =  p  {a [/fP  -  I}'/*’]"  x""' 

+0{K  -  1)  +  P  |7  -  l] 

showing  that  Z  is  again  a  random  variable  with  a  bathtub  hazard  rate 
and  is  easily  generated  as  a  minimum  of  two  independent  Weibull  random 
variables  and  an  independent  exponential  random  variable. 

The  form  given  by  Gaver  and  Acar  (1979)  uses  a  Pareto  distribution  for 
the  early,  decreasing  hazard  rate  instead  of  the  Weibull(p,a)  here.  As 
we  have  seen  in  (iii)  this  will  also  be  suitable  for  a  minification  process. 
Although,  in  both  the  cases  considered  here,  the  process  is  easy  to  generate 
the  correlation  structure  is  difficult  to  determine  analytically. 
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5  Monotonic  Transformation  of  the  Minification  Pro¬ 
cess 


As  noted  in  the  introduction,  one  useful  way  to  generalize  the  minification  pro¬ 
cess  given  by  (  7)  is  to  start  with  the  negative  exponential  marginal  case,  Tavares 
(1980a,  b),  and  take  a  monotonic  transformation  of  each  X^-  Thus  suppose  g 
is  a  monotone  increasing  function;  then  we  define  Yn  =  ff(X„)  and  =  g(Z„) 
for  each  n.  Recall  that  if  Xn  is  negative  exponential  of  mean  1  then  is  also 
negative  exponential  and  of  mean  (K  —  1)~*. 

It  is  straightforward  to  verify  that  the  process  {Vn}  is  stationary  and  Markov 
and  defined  by 


=  min(<,  [«■<?-'  (Vn-i)]  ,9 
=  g  [Kg-^  min  (Vn-i,  W'n)]  • 


(19) 


If  g  is  monotone  decreasing  then  we  must  replace  min  in  (19)  by  max. 

Note  that  by  definition  Xn  is  negative  exponential  of  unit  mean  so  that  if 
Yn  is  to  have  cumulative  hazard  function  Ay(y)  then 


exp(-Ar(y))  =  SK(y)  =  Sx  (3"*(y))  = 


Hence 

S~‘(y)  =  Ay(y).  (20) 

Note  also  that  if  g  is  decreasing  then  -g  is  increasing  so  we  need  consider  only 
increasing  transformations. 

The  bivariate  distribution  of  any  two  observations  may  be  obtained  from 
(  13)  and  is  given  by  the  joint  survivor  function 


S’y„,K„-y(y.*)  = 


Sy  [max|i,y  (g  ‘(y)/^^)}]  5y(y) 


(21) 


SY{9[9~\y)/Ki)) 

Again,  this  bivariate  distribution  is  mixed  with  probability  given  by  (  14),  i.e. 
Py  =  P[g-^{Yn)  =  K^g~^{Yn-}))  distributed  on  the  curve  Kn  =  9  [k^9~^  (^n-j)]. 
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and  the  remaining  (1  -  pj)  distributed  over  ¥„  <  g  (^n->)]  by  means  of 

the  survivor  function  (  21).  When  g  is  monotonic  decreasing  we  replace  survivor 
functions  by  distribution  functions  in  (  21). 


We  can  derive  the  autocorrelation  function  for  {Vn}  from  first  principles  in 
the  same  way  as  for  {Xn}.  However,  we  can  also  relate  the  moments  of  the 
transformed  series  to  those  of  the  original  negative  exponential.  For  example, 

E{Y„Y„-i)  =  E  |fl(X)  |5(XX)5z(X)  +  g(Kz)fz{z)dz^  |  (22) 

which,  when  ^(0)  is  finite,  can  be  simplified  to 

f;(y„y„_i)  =  e  |(;(x)  |<7(o)  +  k  sf[Kz)Sz{z)d^  |  (23) 

where  X  is  negative  exponential  of  unit  mean,  and  ^^(x)  = 


6  Examples  of  the  Transformation  Process 

(i)  y„  =  -  lnX„  so  that  Yn  has  the  extreme  value  or  Gumbel  (1)  distribution 
with  distribution  function  *.  Note  that  g  is  monotonic  decreasing  and 
so  the  process  is  defined  by 

y„  =  max(yn_i,W„)  +  6, 

where  6  =  -  InX.  In  addition,  the  innovation  process  {W^}  is  also  an 
extreme  value  random  variable.  This  is  exactly  the  process  introduced  by 
Tavares  (1977)  and  examined  in  some  hydrological  contexts.  He  was  unable 
to  specify  the  autocorrelation  function  then  but  noted  it  appeared  to  be 
exponential.  An  examination  of  (  22)  suggest  it  would  be  very  difficult  to 
obtain  in  closed  form. 
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(ii)  y„  =  Xn'^,K  >  0,  80  that  Y„  is  Weibull  with  parameter  k.  We  find  that 
the  process  is  defined  by  (  16),  i.e. 

This  is  in  effect  the  Weibull  process  discussed  by  Sim  (1986),  and  in  an 
earlier  example  in  this  paper. 

(iii)  y„  =  >  0,  so  that  y„  is  the  second  type  of  extreme  value 

distribution  (Johnson  and  Kotz,  1970,  Ch.21)  with  distribution  function 
Fy  (y)  =  exp  (-y“*) ,  y  >  0.  The  process  is  defined  by 

y„- A-*max(y„_,,W„). 

but  again  the  autocorrelation  appears  unobtainable  in  closed  form. 

(iv)  y„  =  1  -  so  that  y„  is  now  uniformly  distributed  on  (0,1).  Here,  we 
have  y~^(y)  =  -  ln(l  —  y)  and  the  transformation  is  an  increasing  one,  so 
that  (9)  may  be  used  directly  to  show  that 

y„  =  min(l  -  (1  -  yn-i)^,  1  -  (1  -  W„)^) 

where  IV„  has  survival  function  Sw{v;)  =  (1  -  The  formulation 

may  be  simplified  by  using  the  decreasing  transformation  y„  =  in 

which  case  we  obtain 

=  (24) 

where  {/„  is  uniformly  distributed  on  (0,1).  Both  processes  share  the  same 
autocorrelation  function  and  application  of  (  23)  yields  it  in  the  form 

PxU)  =  ^  1 

(v)  The  uniform  process  of  (iv)  is  important  also  because  we  may  use  it  to 
generate  any  other  marginal  distribution  by  means  of  the  inverse  distribu¬ 
tion  function  transformation.  Thus,  if  Y  is  to  have  distribution  function 
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Fy{X)  we  can  generate  the  process  using  the  transformation 

y„  =  Fy'  (1  -  .  (26) 

Now,  =  —  In5)'(z),  and  we  note  that  this  is  the  cumulative  hazard 

function  of  Y.  Note  that  (  26)  is  a  monotone  increasing  transformation  so 
that  all  the  results  (19)  -  (23)  may  be  applied  directly.  In  particular,  note 
that  if  (  23)  is  applicable  it  may  be  written  in  the  form 

(  r  ftlKX)  , 

£?(ynyn-i)  =  £?|fl(X)l5(0)  +  {Sy(y))‘~'^dy 

where 

aix)  =  Fy\l-e-‘). 

Note  also  that  the  results  (19)  --  (23)  hold  when  {Kn}  is  a  transformed 
version  of  {X„},  given  by  y„  =  y(X„),  given  only  that  {Jfn}  minifi- 
cation  process  satisfying  (7).  We  have  considered  the  case  where  {Xn}  is 
the  negative  exponential  process  but,  in  fact,  {Xn}  may  be  amy  minifica> 
tion  process  and  (19)  -  (23)  still  hold  with  obvious  modifications  to  the 
comments  immediately  after  (23).  In  particular,  it  may  well  be  simpler 
to  choose  to  be  uniform  on  (0,1)  and  then  take  g{x)  =  Fy^{x).  We 
have  specified  two  distinct  uniform  minimum  processes  in  this  paper  and 
either  may  be  used  in  this  way,  although  the  first,  example  (ii)  of  Section 
4,  is  much  simpler  to  implement. 

As  an  example  of  such  a  procedure  we  consider  the  distribution  of  example 
(i)  of  Section  4.  We  noted  that  for  .K’c(l,3)  no  minification  process  exists 
with  this  marginal  distribution  .  We  show  here  how  to  derive  a  suitably 
transformed  process  beginning  with  the  uniform  minification  process  of 
example  (ii)  in  Section  4.  Now, g{x)  =  Fy^{x),  i.e. 

2x  0  <  I  <  0.5 

10*  -  4  0.5  <  *  <  1. 
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For  simplicity,  we  take  K  =  2,  in  which  case  g  \Kg  ^(i)]  =  g{x)  and  so 
the  process  is  given  by 


V„  =  min(y„_i,W^„) 


y  2V„  0<V„<  0.5 

\  lOVn  -  4  0.5  <  <  1, 


where  fVn  =  2Zn  =  2(7„/(l  +  £/«),  from  example  (ii)  of  Section  4,  and  l/n 
is  uniform  on  (0,1). 

7  Maximum  Processes 

We  may  define  a  maximum  process  {Xn}  in  exactly  the  same  way  as  the  mini- 
hcation  process  of  (7).  We  use 

X„  =  amax(X„_i,X„)  0<a<l  (27) 

and  now  the  analog  of  (8)  replaces  survival  functions  by  distribution  functions 
and  an  argument  similar  to  the  one  about  fC  shows  ae(0, 1).  The  discussion 
goes  through  in  an  exactly  analogous  fashion  but  is  somewhat  less  elegant  since 
there  is  no  natural  interpretation  of  the  analogue  of  the  hazard  rate  when  using 
distribution  functions.  Nevertheless,  most  results  are  duplicated  with  survivor 
functions  replaced  by  distribution  function.  For  example,  the  analog  of  (17)  is 
exactly  the  same  with  S^(z)  repleiced  by  F^(z)  =  P(Z  <  z). 

As  a  simple  example  we  present  the  uniform  maximum  process.  Since 
Pz(^)  =  Px(o‘^)/Px(^)  it  follows  that 


_  j  0  with  probability  a 

~  I  a  +  (1  -  a)U  with  probability  1  -  a  ' 

where  U  is  uniform  on  (0,1).  We  may  also  derive  the  autocorrelation  function 

for  this  process  and  it  is  in  the  geometric  form  of  an  AR(1).  In  particular 

Px(i)  =  af*'  i  =  0,l,....  (29) 
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Note,  however,  that  because  of  the  special  mixed  form  of  Zn  we  can  rewrite  the 
process  (  27)  in  the  form 


a  +  (1  -  a)Un 


with  probability  a 
with  probability  1  -  a 


(30) 


Obviously  this  result  is  specific  to  the  uniform  maximum  but  it  serves  to  illus¬ 


trate  the  differences  which  can  arise  between  maximum  and  minimum  processes. 


8  Numerical  Example:  Three  Uniform  Processes 

Three  uniform  autoregressive  process  have  been  derived  in  this  paper.  The  first 

is  the  minification  process  (ii)  of  Section  4,  with  f>x{3)  —  {l/^)^>  i  =  0»  _ 

The  second  uniform  process  is  the  transformation  process  given  as  example  (iv) 
of  Section  6,  with  correlation  given  in  (  25)  as  px(l)  =  3/{2K^  -|-  1).  The  third 
uniform  process  given  is  the  maximum  process  given  in  Section  7  with  Z  given 
at  (28). 

Sample  paths  are  shown  for  these  three  processes  in  Figure  1.  All  three 

are  generated  from  the  same  uniform  i.i.d  sequence  Un,  n  =  0,1 . 100.  Also 

all  three  have  the  same  value  for  Px(l),  namely  0.9.  This  means  that  for  the 
minification  process  K  =  \/px{l)  ~  l/(0.9),  but  for  the  transformation  pro¬ 
cess  K  =  7/6.  For  the  third  process  a  =  (0.9)°®,  from  (29).  Note  also  that 
the  marginal  distributions  are  uniform  and  therefore  bounded  by  zero  and  one, 
unlike  most  time  series  for  which  sample  paths  or  data  are  exhibited. 

In  the  top  panel  of  Figure  1,  we  see  the  typical  “runs  up”  or  “run  off”  behavior 
of  a  minification  process.  Again  the  middle  panel  shows  that  the  transformation 
process  exhibits  “runs  down” ,  with  the  runs  being  convex  down  when  they  start 
at  high  values,  and  convex  up  when  they  start  at  low  values.  Note  particularly 
the  very  slow  decay  at  the  end  of  the  series  when  the  process  has  very  high  values. 
In  the  third  panel,  the  maximum  uniform  process  shows  very  long  geometric 
decays  and  has  a  very  odd,  persistent  appearance. 
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Figure  2  shows  the  autocorrelation  functions  of  the  transformation  and  the 
minihcation  uniform  processes.  The  function  px{j)  minification  process 

is  less  than  the  function  pxU)  transformation  process.  In  fact,  for  large 

jf,  px{j)  for  the  transformation  process  has  value  approximately  one  and  half 
times  that  of  the  autocorrelation  pxU)  for  the  minification  process. 
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Figure  1:  Sample  paths  for  three  processes  with  uniform  marginal  distributions. 
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Autocorrelation 


Uniform  Processes 

Autocorrelotion  Functions  for  two  Processes 


Figure  2:  Correlation  functiona  Px(j)  for  the  minification  and  transformation 
processes  with  uniform  distribution. 
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